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Abstract 

We predict that the four-loop contribution 03 to the QCD (3 function in the MS 
prescription is given by (3 3 ~ 23,600(900) - 6, 400(200) iV) + 350(70) Nj + 1.5 Nf, 
where Nf is the number of flavours and the coefficient of Nj is an exact result from 
large- Nf expansion. In the phenomenologically-interesting case Nf = 3, we estimate 
(3 3 = (7.6 ± 0.1) x 10 3 . We discuss our estimates of the errors in these QCD predictions, 
basing them on the demonstrated accuracy of our method in test applications to the 
O(N) $ theory, and on variations in the details of our estimation method, which goes 
beyond conventional Pade approximants by estimating and correcting for subasymptotic 
deviations from exact results. 
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A perennial problem in Quantum Field Theory is the search for calculational 
techniques that go beyond conventional perturbation theory. These are needed in 
order to control higher-order terms in what are in many cases, such as QED and 
QCD, asymptotic perturbation series, and to make contact with essentially non- 
perturbative methods such as the lattice. To be credible, any proposed technique for 
going beyond perturbation theory must demonstrate its ability 1) to make significant 
predictions with 2) reliable error estimates. One candidate technique is that of Pade 
approximants, which have previously met these criteria in applications to problems 
in condensed- matter theory and statistical mechanics 0, in particular. 

We have been working to extend these successes to four-dimensional Quantum 
Field Theories such as QED and QCD However, the perturbative calcula- 

tions at any given order in these gauge theories are considerably more complicated 
than in many lower- dimensional field theories. Thus the calculated perturbative 
QCD series are often much shorter than those familiar in condensed- matter and 
statistical-mechanical applications. This is on the one hand an opportunity, be- 
cause it accentuates the need for auxiliary techniques. However, it is on the other 
hand a challenge for the Pade approach, which is based on the available perturba- 
tive results. Shorter series provide less input into the Pade machinery, increasing 
unavoidably the errors and making it difficult to cross check their magnitudes. 

A welcome opportunity to test the Pade method is currently provided by the 
perturbative series for the renormalization-group (3 function in QCD Q, which is 
particularly interesting for several reasons: 

1) It is the most fundamental quantity in perturbative QCD, and its behaviour 
impinges on many of the most basic issues in this theory, such as the question 
whether the QCD coupling strength a s approaches a finite value in the infrared 
limit §. 

2) It has already been calculated exactly to three-loop order for any number 
of flavours Nf j|, providing as much information as is available from any 
perturbative QCD series to be incorporated into any attempted extrapolation. 

3) This and the availability of four- and five-loop calculations of the (3 function 
in O(N) scalar field theory in four dimensions H,[[]|§, which has a similar 
structure, facilitate the estimation of plausible errors. 

4) The calculation of the four-loop contribution (3% to (3 seems to be within the 
scope of available perturbative techniques, and we are aware of at least one 
project || that is under way to evaluate it, which means that a prediction can 
be checked in the foreseeable future. 

We predict in this paper that 

& ~ 23, 600(900) - 6, 400(200) N f + 350(70) Nj + 1.5 iVj (1) 

in the MS prescription, where the coefficient of the last term is an exact input cal- 
culated using large- Nf techniques [fLOR . In the case Nf = 3 of phenomenological 
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interest, we predict (3$ = (7.6 ± 0.1) x 10 3 . The estimates of the errors in our pre- 
dictions are discussed in more detail below. They are based on the demonstrated 
accuracy of our method in test applications to the O(N) $ 4 theory, and on varia- 
tions in the details of our estimation method, which goes beyond conventional Pade 
approximants by estimating and correcting for subasymptotic deviations from exact 



results [ITT 



We now review the formalism of Pade approximants |f[2 |, introducing the nota- 

oo 

tion we use in this paper. Considering a general perturbative series P(x) = ^ S n x n , 



n=0 



the Pade approximant P[n/m](x) of order [N/M] is the ratio of polynomials A^{x) 
of order N and Bm(x) of order M, chosen so that 



[N/M] 



A N {x) 
B M (x) 



P(x) + 0(x 



N+M+l^ 



(2) 



One may use the P\n/m]{x) in either of two ways: one may interpret the coefficient 
£]v+m+i °f x N+M+1 in a power-series expansion of P[n/m]( x ) as the Pade Approx- 
imant Prediction (PAP) for the next perturbative coefficient Sjv+m+Ij or one may 
interpret the full expression P\n/m}{ 

X^) ctS cl Pade Summation (PS) of the entire series 
P(x) ||. In this paper, we shall mainly be concerned with refinements of the PAP's 
for Quantum Field-Theoretical perturbation series, but we shall also comment at 
the end on an application of PS. 

We have exhibited previously 0[| mathematical conditions sufficient to guaran- 
tee that PAP's converge to Sn+m+i as N and M are increased. This is in particular 
true for any series whose perturbative coefficients have the property that 



Sn+lSn-l 



or 



A 

1 + - 

n 



(3) 



as n — > oo. The first of these conditions is satisfied by any perturbative series that 
is dominated by a finite set of renormalon poles, and the latter condition is satisfied 
by the Borel transforms of such series, guaranteeing the utility of PAP's in many 
applications to perturbative QCD [El. 



We have, moreover, been able to derive 
sponding relative errors of PAP's: 



asymptotic expressions for the corre- 
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N+M+l 



oPade c 
o N+M+1 — D N+M+l 



s 



N+M+l 



M\A M 
N M 

M\B(B + l)...(B + M - 1) 

N 2M 



(4) 
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as iV — > oo, for fixed M. There is empirical evidence that the relative errors of 
PAP's are fitted even better by expressions of the form (H) with the replacements 

N^N + M + aM + b (5) 

where the parameters a, b are to be fitted on a case-by-case basis.0 See for example 
the comparison between PAP's and the true values of perturbative coefficients in 
QCD for a large number of flavours Nf — > oo shown in and compared with 



asymptotic formulae for the relative errors in [13|] 



These convergence theorems and the successes of the asymptotic error formulae 
(fD lead us to propose that one use the estimated asymptotic errors to correct PAP's 
systematically to the Asymptotic-Pade Approximant Predictions (APAP's): 

cPade 

qAPAP °N+M+1 /n\ 

^N+M+l 1 i A W 
1 + On+M+1 



We document elsewhere [ I I] numerical evidence that APAP's excel even over con- 
ventional PAP's in the accuracy with which they reproduce known coefficients in 
QCD perturbation series, using as inputs calculations of lower-order coefficients. 
Our objective in this paper is to predict, with a meaningful and credible error es- 
timate, the next unknown (four-loop) term in the perturbative series for the QCD 
(3 function, whose exact calculation by conventional methods is currently being un- 
dertaken §|. If our prediction turns out to be correct within the stated errors, it 
will provide reason to believe Pade predictions for other coefficients that are less 
tractable to calculate exactly. 

We denote the QCD (3 function by f3(x) = - (J2n=o (3 n x n+1 ), where x = a s /(An). 
The known coefficients for N c = 3 in the MS prescription are || (3q = 11 — (2iV//3), 
/3i = 102 - (38N f /3) and f3 2 = (2857/2) - (5033/18)]^ + (325/54)iV|. It is pos- 
sible to use the [0/1] PAP to estimate /3 2 on the basis of the exact values of /? ,i- 
This is done by fitting the numerical values of /^{Nf) obtained from the [0/1] PAP 
for Nf = 0, ... ,4 to the polynomial form /?2 = A + BNf + CNj. Since the leading 



coefficient of Nf is known at each order in loop expansion |]T0| , we anticipate here 
the approach to be used in the 4-loop case and include its exact value in the fit, 
finding 

(3 3 = 946 - 177 N f + 6.02 N 2 f (7) 

Although qualitatively correct, this estimate is not very accurate, which was hardly 
to be expected. The use we make of the [0/1] PAP is to estimate, for each individual 
value of Nf, the normalization coefficient A or B of the asymptotic correction (U), 
which we then use as an input in calculating the [1/1] APAP @. 

The result we obtain clearly depends on our choice between the possible asymp- 
totic forms in @, which is /3 n+ i/3 n _i//^ ~ 1 + A/n. In contrast to the case with 

Jin certain simple cases such amended formulae for the relative errors (0) turn out to be exact. 
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many QCD observables such as an/ir = [-R e+e _ — J2fQ}) i n e+ e ~ ~~ > hadrons, or 
the anomalous dimension of the quark mass insertion, we are unaware of any proof 
that there is a renormalon: ft n ~ n\ in the MS prescription used in this paper, which 
would be a sufficient condition for this assumption to apply || g. More is known 
about the ft function in supersymmetric QCD. Indeed, it is trivial in N = 2 QCD, 
being given entirely by the one-loop term |TJ|. In the case of N — 1 QCD with 
matter, it is known that in a background-field Pauli-Villars prescription [15[ 



8(a) = 9 " ?>N c -N s + 1 { g 2 )N } 

^ Ky> 16tt 2 l-N c (g 2 /8ir 2 ) K> 

where 7(g 2 ) is the anomalous dimension of the mass operator. If the perturbative 
coefficients in 7(g 2 ) are dominated by a renormalon, then also ft n ~ n\, and our 
asymptotic assumption applies []. On the other hand, if one considers the limit 
N c g 2 fixed and Nf/N c = 3 — e fixed, the numerator of /3(g) fl8|) is determined by the 
0(g 2 ) term in 7, and the [1/1] PAP is exact] 

These considerations support our assumption that the leading behaviours of 
the perturbative coefficients are as given in the top line of (§) , leading to the corre- 
sponding expression (^) for the relative error, and hence justifying the corresponding 
APAP correction @, though they by no means prove it. 

As a cross check on our method, and to aid in assessing its likely errors, we first 
apply it to the O(N) $ 4 theory in 4 dimensions, whose ft function we denote by 
ft (9) — J2^=o ftng n+2 '■ These ft function coefficients are known for n < 4, and it is 
known that ft n ~ n\ as n — > 00 in this theory ]TE |, justifying the APAP procedure 



with the first form of subasymptotic correction in (|4]). We follow the same procedure 
as advocated above for QCD, namely we use the known results for /3 0j i in [0/1] PAP's 
for various values of N and then apply the APAP correction (|6]) to the naive [1/1] 
PAP |. Our results for ft 3 in O(N) $ 4 theory are shown in Fig. 1(a). We see that 
the naive [1/1] PAP falls considerably below the known exact result for ft 3 in this 
case, which is known to be a polynomial in iV: (3% = A + BN + CN 2 + T>N 3 , shown 
as the solid line in Fig. 1(a). It is clear from this figure that the [1/1] APAP with 
the A/n type of correction (||) is startlingly accurate, reproducing qualitatively the 
known polynomial form and quantitatively the values of ft 3 for N < A. We have 
in fact determined empirically that the most accurate prediction of the coefficients 
in the polynomial expansion of ft% may be obtained by a naive average (A) of the 
correction coefficient A over the values iV = 0,1,2,3,4 studied. Using this (A) 

^Consider, however, a generic change in renormalization prescription from MS to one in which 
the QCD correction to an observable (e.g., or) is treated as an effective charge. In order to 
cancel the factorial growth of the MS perturbative coefficients for the observable, it is clear that 
the /3-function coefficients must also diverge factorially. 

^We note that the positions of renormalon singularities are not thought to depend on the 
renormalization scheme chosen, just their strengths. Therefore, one may expect this conclusion to 
carry over to the MS prescription used in this paper. 

" This was used as an order-of-magnitude estimate in ]l7[ . 
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prescription, and then making a polynomial fit to (3% = A + BN + CN 2 + DN 3 , with 
the exactly known value of the 0(N 3 ) coefficient T> ||[7|||, we find the numerical 
values of the polynomial coefficients shown in Table I, reported together with their 
errors relative to the known exact coefficients For comparison, we also display 
the relative errors found if one fits the correction coefficient A individually for each 
value of Nf, as well as the errors found making the (disfavoured) B/n 2 correction, 
and the errors of the naive PAP. 



Table I 

APAP results for the polynomial coefficients of (3% in O(N) $ 4 



polynomial 


exact 


APAP 


rel. error 


rel. error 


rel. err. 


rel. err. 


coefficient 


calculation 


prediction 


APAP 


APAP, 


APAP, 


PAP 








(A) 


A(N f ) 


(B) 




A 


-100.46 


-104.32 


+3.84% 


+4.69% 


-17% 


-31% 


B 


-33.28 


-34.39 


+3.34% 


+1.38% 


-17% 


-31% 


C 


-2.06 


-1.67 


-18.87% 


-17.18% 


-35% 


-45% 


V 


6.4 x 10~ 4 


input 











We note that the APAP's with the B/n 2 form of correction, also shown in 
Fig. 1(a), are considerably less accurate than those with the A/n form, in accord 
with what we expect from the asymptotic results of |TB| . We also note that the naive 



[0/1] PAP is a significantly worse approximation to 02 in O(N) $ 4 theory than the 
corresponding [0/1] PAP in QCD. This implies that the APAP correction is rela- 
tively larger in this case, and leads us to suspect that the [1/1] APAP prediction in 
QCD might even be more accurate than that in Fig. 1(a) or Table I. 

We now make our prediction for the 4- loop (3 function coefficient (3$ in QCD. As 
in the previous O(N) <3> 4 case, we first calculate the [1/1] PAP's for < Nf < 4, 
and then make the appropriate APAP correction assuming that (j3 n+ ij3 n -i) / /3 2 ~ 
1 + A/n, obtaining the predictions shown as crosses in Fig. 1(b). These points are 
obtained by averaging the value of the asymptotic coefficient A over the values of 
Nf chosen, which is known to give accurate predictions in the O(N) $ 4 case, as 
we have already commented in connection with Fig. 1(a). In the particular case 
Nf = 3 of phenomenological interest, we predict (3% = (7.6 + 0.1) x 10 3 , with an 
error which we estimate conservatively from the corresponding value of N in the 
$ 4 model, and for N f = 1,2,4 we find (3 3 = (17.6,12.3,3.8) x 10 3 with similar 
fractional errors. We also fit the APAP's to the expected polynomial dependence of 
form (3 3 = A + BN f + CNf + VNf, where the O(Nf) coefficient V is known exactly 
|jiT||, obtaining the predictions shown in Table II. 



**The errors in the polynomial coefficients are correlated in the fitting procedure, and our errors 
for individual values of iV are smaller. In the particular case N = 3, for example, we find = 222.4, 
to be compared with the exact result (3% = 218.8, an error of 1.6%. 
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Table II 



APAP results for the polynomial coefficients of (3% in QCD 



coefficient 


APAP, (A) 
prediction 


error 
estimate 


APAP, A(Nf) 
prediction 


PAP 
prediction 


A 


23,557 


±900 


24, 086 


20,015 


B 


-6,353 


±200 


-6, 572 


-5,396 


C 


346 


±70 


357 


292 


V 


1.5 (input) 









We believe on the basis of our 0(N) $ 4 experience that the most reliable pre- 
dictions are provided by the (A) APAP's in Fig. 1(b), with the relative errors of the 
corresponding O(N) $ 4 coefficient predictions providing a conservative estimate of 
the errors in the QCD case. These errors, which we expect to be highly correlated 
as in the $ also shown in Table II, as well as the APAP's obtained using 

A(Nf). We prefer the (A) APAP's, but are pleased to note that the differences 
are comparable with our quoted error estimates. As an exercise, we have also com- 
puted and shown in Fig. 1(b) as open circles the APAP results obtained using the 
B/n 2 form of correction, and also display the naive PAP's shown as diamonds. We 
believe that the differences from the B/n 2 APAP's are likely to be overestimates 
of the errors: as we have discussed above, the B/n 2 APAP's have less theoretical 
motivation. Finally, Table II also displays the numerical differences between our 
preferred APAP's and the naive PAP's. 

As already mentioned, we are aware of a project underway to calculate /3 3 exactly 
in QCD, which will enable our predictions to be tested. If the exact results are 
consistent with the predictions of our APAP method, we will then be able to use 
them to predict the 5- loop coefficient /3 4 in QCD. Fig. 2 shows the corresponding 
(A) and (B) predictions for /3 4 in the O(N) $ 4 model, compared with the known 
exact result in this model [J. The corresponding (A) predictions for the polynomial 
coefficients of /3 4 for this model are displayed in Table III. They were obtained by 
fitting (A) and the combination a ± b of subasymptotic coefficients in (^) to the 
known exact results for {3% and (3$ in the O(N) <3> 4 model. Also shown in Table III 
are the known exact results and the relative errors. We see in Fig. 2 that the (A) 
and (B) predictions are closer in this higher-order case, as one might have expected. 
For comparison, we recall the prediction |7j that /3 4 = 1405 ± 80, which is to be 
compared with the exact result /3 4 = 1424.29 || and our prediction that /3 4 = 1432, 

^We note in passing that our APAP approach easily revealed the sign misprint that appeared in 
the original version of , and, even after sign correction, agrees better with the corrected value || 
than with the previous value j?J that had a small numerical error. We recall that PAP's have on 
previous occasions shown that they can identify errors and misprints in tabulations H . 
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corresponding to an error of 0.56%. We would hope that predictions of comparable 
accuracy will be possible in QCD, once the exact result for (3% is available. For 
completeness, we also provide our APAP prediction for the 6-loop coefficient f3$ in 
the O(N) model, which has not yet been computed exactly: for N = 0,1, 2, 3, 4, we 
predict -fc = 11828, 17687, 24958, 33802, 44330, respectively. For comparison, we 
recall the prediction that -f3 5 = 17, 200 ± 50 in $ 4 theory with N = 1. 



Table III 

APAP results for the polynomial coefficients of (3± in O(N) $ 4 



coefficient 


exact 


APAP 
prediction 


rel. error 


A 


1002.0 


1002.26 


0.03% 


B 


385.6 


390.61 


1.3% 


C 


36.12 


38.82 


7.5% 


V 


0.576 


0.627 


8.8% 


£ 


-0.0013 


input 





It is not the purpose of this paper to discuss in detail the possible phenomeno- 
logical implications of our prediction for f}$ in QCD []18| . We limit ourselves here 



to advertizing the fact that the conventional PS procedure, applied to our APAP 
result for Nf, predicts the existence of a zero in the QCD (3 function, and hence an 
infrared fixed point in the MS coupling, whose location coincides with that predicted 
previously by Stevenson || using very different arguments. 
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Note added in proof: 

After this paper was submitted, T. van Ritbergen, J. A.M. Vermaseren, and 
S.A. Larin presented the results of an exact calculation of the four-loop QCD j3 
function ( |hep-ph / 97013901) . The result in their eq. (8) contains qualitatively new 
color factors, corresponding to quartic Casimirs, which are not present at 1-, 2- 
and 3-loop level, and therefore cannot be estimated using the Pade method. Con- 
sequently, our predictions should be compared with the rest of their exact ex- 
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pression, as shown below. Writing (3% = A + BNp + CNp + VNp, we have: 



coefficient 


APAP estimate 


exact 


relative 


number of 




(uncertainty a) 


result 


error 


a's 


A 


23,557(900) 


24,633 


-4.4% 


1.20 


B 


-6,353(200) 


-6375 


-0.35% 


0.11 


C 


346(70) 


398.5 


-13.2% 


0.75 


V 


1.5 


1.5 


input 
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Figure 1: (a) The 4- loop (3- function coefficient in $ 4 theory with O(N) symmetry. 
The exact results are denoted by black dots, joined by a solid line to guide the eye. Naive 
PAP results are denoted by diamonds, and APAP results obtained from the A/n type of 
correction are denoted by crosses. For comparison, also shown are APAP results obtained 
from the B/n 2 type of correction, denoted by open circles, (b) The 4- loop (3- function 
coefficient (3% in QCD with Nf flavours, using the same notation as in Fig. 1(a). The A/n 
APAP results are assigned error bars obtained from the relative errors in the corresponding 
<£ 4 results. 
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